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ABSTRACT 

We discuss an analytic approach for modeling structure formation in sheets, 
filaments and knots. This is accomplished by combining models of triaxial col- 
;h ■ lapse with the excursion set approach: sheets are defined as objects which have 

collapsed along only one axis, filaments have collapsed along two axes, and halos 
are objects in which triaxial collapse is complete. In the simplest version of this 
approach, which we develop here, large scale structure shows a clear hierarchy of 
morphologies: the mass in large-scale sheets is partitioned up among lower mass 
filaments, which themselves are made-up of still lower mass halos. Our approach 
provides analytic estimates of the mass fraction in sheets, filaments and halos, 
and its evolution, for any background cosmological model and any initial fiuctua- 
tion spectrum. In the currently popular ACDM model, our analysis suggests that 
more than 99% of the cosmic mass is in sheets, and 72% in filaments, with mass 
larger than 10^°Mq at the present time. For halos, this number is only 46%. Our 
approach also provides analytic estimates of how halo abundances at any given 
time correlate with the morphology of the surrounding large-scale structure, and 
how halo evolution correlates with the morphology of large scale structure. 



Subject headings: structure formation: ellipsoidal collapse - mass functions: 
sheets, filaments, halos 

1. Introduction 

Recent observations and high resolution numerical simulations of structure formation show 
that, on large scales, the Universe is best thought of as a Cosmic Web (Bond, Kofman & 
Pogosyan 1996): large scale sheets are traversed by filaments, which themselves intersect at 
knots (Bhavsar & Ling 1988; Bharadwaj et al. 2000; Bhavsar et al. 2003; Sheth et al. 2003; 
Bharadwaj & Biswajit 2004; Shandarin ct al. 2004; Biswajit & Bharadwaj 2005; Colberg et 
al. 2005), and there are vast relatively empty voids in between (Gregory & Thompson 1982; 
Kirshner et al. 1983; Miiller et al. 2003; Vogeley et al. 2004). The knots where filaments 
intersect are the locations of rich galaxy clusters, and these are often identified with the 
massive virialized dark matter halos found in simulations. Halos are the best-studied features 
of the cosmic web; they tend to be about two hundred times denser than the background 
universe, so they account for a small fraction of the volume of the Universe. Simple analytical 
models have been very successful in understanding the properties of dark halos (Gunn & Gott 
1972; Press & Schechter 1974; Bond et al. 1991). However, with the exception of Zeldovich 
(1970) and Gurbatov, Saichev & Shandarin (1989), a framework for describing the properties 
and dynamics of filamentary and sheet-like structures in the cosmic web is still lacking. The 
main goal of this paper is to provide such a framework: our goal is to discuss the cosmic web 
using language which extends naturally that currently used for discussing dark halos. This 
complements recent work showing how the same framework can be used to discuss cosmic 
voids (Sheth & van de Weygaert 2004). 

In current descriptions of large scale structure, the abundance by mass of dark matter halos 
and the evolution of this abundance, i.e., the cosmological mass function n{M, z), plays a fun- 
damental role. Models of this mass function suggest that it encodes information about both 
gravitational dynamics and the statistics of the initial fluctuation field (Press & Schechter 
1974; Bond et al. 1991). Models of n{M,z) which assume that dark halos form from a 
spherical collapse are in reasonable but not perfect agreement with results from numerical 
simulations of hierarchical gravitational clustering (Sheth & Tormen 1999). In reality, halos 
are triaxial (Jing & Suto 2002), so they cannot have formed from a spherical collapse. The 
initial shear field and tidal effects almost certainly play some role in determining the evolu- 
tion of an object, and models of the associated non-spherical collapse have been developed 
(Zeldovich 1970; Icke 1973; White & Silk 1979; Peebles 1980; Bond & Myers 1996). Such 
non-spherical collapse models can be incorporated into models of the halo mass function in 
various ways (Eisenstein & Loeb 1995; Monaco 1995; Lee & Shandarin 1998; Chiueh & Lee 
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2001). In particular, Sheth et al. (2001) argued that much of the discrepancy between model 
predictions and simulations is removed if one replaces the assumption of a spherical collapse 
with one where halos form from a triaxial collapse. 

In this paper, we exploit a crucial difference between spherical and ellipsoidal collapse models 
of halo formation. Namely, if we think of a halo as being triaxial, then the formation of an 
object corresponds to the time when gravity has caused all three axes to turn around from 
the universal expansion and collapse (Sheth et al. 2001). If the object is spherical, this 
collapse occurred at the same time for all three axes. In contrast, a triaxial object has three 
critical times, corresponding to the collapse along each of the three axes — the shortest axis 
collapses first, the intermediate axis later and the longest axis last. 

If we identify 'halos' with objects which have collapsed along all three axes, then it is natural 
to identify filaments with objects which have collapsed along only two axes, and sheets with 
objects which have only collapsed along one axis. Given this identification, the remainder 
of this paper describes how to construct a model which naturally encapsulates the idea that 
halos form within filaments which themselves populate sheets. In particular, we provide a 
framework for discussing the mass functions of sheets and filaments, and how these evolve, 
as well as the conditional mass functions of halos within filaments and sheets, and their 
evolution. 

Section 2 presents our model. It begins with a discussion of our triaxial collapse model 
(Section 2.1). Appendix A compares this model with others in the literature. It then 
translates this model into a form which is most convenient for use in the 'excursion set' 
approach (Sections 2.2 — 2.4). The excursion set approach is commonly used to estimate 
the mass function of halos (Bond et al. 1991; Sheth et al. 2001). Our collapse model, 
when combined with the excursion set approach, has a rich structure. We illustrate this by 
using it to generate analytic approximations to the mass functions of sheets and filaments 
(Section 2.5), and to the conditional distributions of filaments in sheets, and halos in filaments 
(Section 2.6). Comments on how mass accretion may depend on environment (Section 2.7), 
and on the characteristic densities of sheets, filaments and halos, which may aid identification 
of such objects in the cosmic web (Section 2.8) conclude this section. Section 3 summarizes 
our findings. 
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2. A model for sheets, filaments and halos 

2.1. Ellipsoidal collapse 

The gravitational collapse of homogeneous ellipsoids was studied by Icke (1973), White & 
Silk (1979), and Peebles (1980). Although the exact evolution must be solved-for numerically. 
White & Silk (1979) provide an elegant analytic approximation for the evolution which is 
remarkably accurate. Unfortunately, these early analysis did not reduce, as they should, to 
the Zeldovich approximation in the linear regime. Bond & Myers (1996) noted that this was 
because they did not include the effects of the external tide self-consistently; once these effects 
are included, the collapse model is indeed self-consistent. Appendix A provides details, and 
also shows how the analysis of White & Silk (1979) can be extended to this self-consistent 
case. In this model, the initial collapse is expected to be dominated by the local strain 
tensor, which includes both internal and external tidal forces. For any given cosmology, the 
evolution of a patch in this model is determined by three quantities: its initial overdensity, 
5, and two shape parameters, e and p. It is this model which we will use in what follows. 



2.2. Mass functions and random walks: Three barriers for three cLxes 

In the 'excursion set' approach (Bond et al. 1991), an approximation for the mass fraction in 
bound virialized halos of mass m is obtained by mapping the problem to one which involves 
the first crossing distribution of a suitably chosen barrier by Brownian motion random walks. 
The choice of barrier is set by the collapse model and by the epoch for which one wishes 
to estimate the mass function, and the mapping between random walk variables and halo 
masses is set by the shape of the power spectrum of the initial fluctuation field. 

In particular, if /{u) denotes the fraction of walks which first cross a barrier at scale u, where 
u is the random walk variable, then the mass fraction in objects of mass m is 



„, , 2^("^>-^) d\nm 

(Bond et al. 1991). Here (T^(m) is the variance in the initial density fiuctuation field, 
smoothed on scale R = {3m/4:npmoy^^ , and extrapolated using linear theory to the present 
time. 

If the collapse is spherical, then the barrier is particularly simple: it has constant height, 
with higher barriers required to model the halo population at higher redshift. In particular, 
this height is very simply related to the initial overdensity required for spherical collapse 
by redshift z. This value, which is usually denoted Ssc{z), will play a special role in what 
follows. 
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In the ellipsoidal collapse model, the barrier height associated with halo formation by redshift 
z depends on three numbers: the initial overdensity, and the shape parameters e and j). In 
principle, then, the problem of estimating the mass function is one of crossing a barrier 
in a higher dimensional space (Sheth et al. 2001; Chiueh & Lee 2001; Sheth & Tormen 
2002). Sheth et al. (2001) suggested that the computational complexity could be reduced 
significantly if one used suitably chosen representative values of e and p. Once these have 
been set (for reasons given in Appendix A, they suggested p = and e = [(T(m)/(5]/\/5), 
the excursion set approach corresponds to finding the first crossing distribution by a one- 
dimensional random walk in 5 of a barrier whose height is mass- dependent. In this respect, 
the estimation of the mass function associated with eUipsoidal collapse proceeds exactly as 
for spherical collapse — the only complication is that the barrier height is no longer constant. 

Our model for the mass function in filaments and sheets is entirely analogous. In principle, 
filaments at z can form from some combination of S, e and p; the precise combination depends 
on the details of the collapse model. In our preferred model, the same one used by Sheth 
et al. (2001), this requires that collapse occurs along the intermediate length axis at z, but 
along the longest axis at a later time. Thus, estimation of the mass fraction in filaments 
requires solution of a high-dimensional random walk, and, as for collapse along all three axes, 
we approximate by setting e and p equal to representative values, and considering a one- 
dimensional walk in the remaining variable S. In this case, the barrier height associated with 
filament formation is that associated with collapse along two, rather than three axes. And a 
similar argument means that the mass function in sheets can be approximated by studying 
the first crossing distribution of the one-dimensional barrier associated with collapse along 
just one axis. 

Figure 1 illustrates the basic idea. The jagged fine shows a random walk, and the three 
dotted fines show the barriers associated with collapse along one (bottom), two (middle) 

and three (top) barriers. In our model, the fraction of walks which first cross the lowest 
barrier on scale a{m) equals the mass fraction in sheets of mass m, the fraction which first 
cross the second barrier at cr{m) equals the mass fraction in filaments of mass m, and the 
fraction which first cross the highest barrier at a{m) equals the mass fraction in triaxial 
halos of mass m. We will provide estimates of these quantities shortly. 

However, Figure 1 shows that this model can also be used to provide significantly more 
information. For instance, notice that all walks first-cross the barrier associated with sheets 
at a higher mass scale than when they first-cross the barriers associated with filaments 
or halos. Thus, in addition to providing a way to estimate the mass functions of sheets, 
filaments and halos, our approach also provides a framework for discussing the mass fraction 
in halos of mass that are embedded in filaments of mass m/ < which are themselves 
surrounded by sheets of mass > ruf. Figure 1 shows that this is simply related to the 
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Fig. 1. — An example of a random walk (solid line) crossing the barriers (dotted lines) 
associated with sheets, filaments and halos (bottom to top). The fraction of walks which 
first cross the lowest barrier at ij{ms), then first cross the second barrier at (riiJif) and finally 
cross the highest barrier at a{mh) represents the mass fraction in halos of mass nih which 
are embedded in filaments of mass rrif > rrih, which themselves populate sheets of mass 
TTis > (recall that cr is a decreasing function of m). The precise barrier shapes depend 
on the collapse model; the dotted curves show the barriers in equation (2). 

fraction of walks which are conditioned to pass through certain points before first crossing 
ruh- Thus, our approach allows one to estimate how halo abundances correlate with the 
morphology and density of the large-scale structure which surrounds them. 

So far, we have only discussed the model at a fixed redshift. Since the barrier shapes depend 
on redshift, our approach also allows one to estimate how sheet, filament and halo abundances 
evolve. Moreover, it allows one to address questions like: is the halo abundance in sheets 
of mass lO^^M© at 2; = 1 any different from the mix of halos in lO^^M© filaments today? 
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Our model naturally incorporates the fact that the mass in high redshift sheets becomes 
partitioned up among filaments at later times, with filaments themselves being partitioned 
into halos. Clearly, the model is rich in spatial, temporal and spatio-temporal information. 



The key output from the triaxial collapse models is an estimate of the typical overdensity 
required for collapse along one, two and three axes by redshift z. The dotted curves in 
Figure 1 show how these three 'barriers' depend on mass. Prom bottom to top, the curves 
show 



These analytic approximations to the barriers associated with collapse along one, two and 
three axes are reasonably accurate (c.f. Figure 7). They show that the critical overdensity 
for ellipsoidal collapse along all three axes, (5ec3) is larger than in the spherical collapse model 
(for which this number is ^^c)- However, collapse along just two axes requires an overdensity 
(5ec2 which is almost exactly that in the spherical model, and collapse along one axis only 
requires a smaller initial overdensity, 5eci. In this model, tidal forces enhance collapse along 
the first axis and delay collapse along the last axis relative to the spherical collapse model 
(Sheth et al. 2001) — the expressions above quantify these effects. Note that the differences 
among the three critical overdensities are larger for larger values of cr, corresponding to 
ellipsoids of lower masses. 

Notice that these barrier shapes depend both on aim) and on bgciz). The presence of 
these two terms reflects the fact that the collapse depends on the expansion history of the 
universe, and on the initial spectrum of fluctuations. What is of particular importance in 
what follows is that the barrier shapes actually depend, not on and o individually, but on 
the scaling variable v = [6sc{z)/(7{m)]'^. In the excursion set approach, this implies that the 
mass functions of sheets, filaments and halos at any given time, in any given cosmology, and 
for any given initial fluctuation spectrum, can all be scaled to universal functional forms. 
Our next step is to provide analytic approximations to these forms. 



2.3. Mass dependence of the three barriers 
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2.4. First crossing distributions 



In the simplest excursion set model, object abundances are related to the first crossing dis- 
tributions by one- dimensional random walks of the moving barriers in equation (2). These 
distributions are easily obtained by Monte-Carlo simulation. (Alternative numerical meth- 
ods are also available, e.g., Zhang & Hui 2006.) The Monte-Carlo results are simplified 
considerably by the fact that the barriers can be expressed in terms of the scaling variable 
— [<^sc(-2^)/c("^)]^- The histograms in Figure 2 show the results. The abundances of sheets, 
filaments and halos are related to these first crossing distributions by equation (1). 
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Fig. 2. — Histograms show the first crossing distributions of the barriers associated with 
collapse along one, two and three axes (histograms which extend to largest v are those for 
collapse along only one axis). Curves show analytic approximations to these distributions 
following Sheth & Tormen (2002). 

The smooth curves show analytic approximations to these first crossing distributions, com- 
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puted following Sheth & Tormen (2002). Namely, when the barrier is of the form 



(3) 



as is the case for all three barriers in equation (2), then the first crossing distribution is 
approximately 




-i^[l+^i/-«]2/2 



27r 



1 + 



l-a + 



a{a — 1) 
2! 



+ 



(4) 



Alternative analytic estimates could be computed by noting that when a — 1/2, then the 
first crossing distribution can be written analytically in terms of sums of parabolic cylinder 
functions (Breiman 1967). Given that these would only provide approximations to the first 
crossing distributions we require (because a ^ 1/2), we have chosen the considerably simpler 
approximation given in equation (4). Figure 2 indicates that these approximate solutions 
are sufficiently accurate to allow analytic estimates of a number of interesting quantities. 



2.5. Abundance by mass and redshift 

Upon assumption of a cosmology [which determines 5sc(-2)] and a fluctuation spectrum [which 
determines a{m)], the flrst crossing distributions shown in Figure 2 can be converted to 
mass functions using equation (1). In what follows, wc assume a spatially flat model with 
cosmological constant Aq = 1 — Oq, where {Vlo, h, as) = (0.3, 0.7, 0.9). 

The dashed lines in the bottom-left panel of Figure 3 show the evolution of the number 
densities of halos having mass lO^^M© (top) to lO^^M© (bottom). Except for the lowest 
mass range, these curves decrease monotonically with increasing redshift. This reflects the 
well-known hierarchical nature of structure formation: at any given redshift, massive halos 
are less common than halos of lower mass, and they are even less common at earlier times. As 
a result, if one integrates the mass fraction in objects above some minimum mass, then this 
fraction is smaller at higher redshift. If, instead, one flxes the mass fraction, then one must 
integrate down to smaller mass halos to obtain this fraction. Since less massive halos are 
more abundant, plots of the abundance associated with a flxed mass fraction increase with 
increasing redshift (Mo & White 2002). The solid curves show such locii in the abundance- 
redshift plane for a range of mass fractions. The top-left and right panels show analogous 
results, but for sheets and filaments respectively. Similar trends are seen in all three panels, 
indicating that sheets and filaments grow hierarchically in much the same way as halos. 

These panels show that, at z = 0, more than half of the cosmic mass is in sheets with masses 
exceeding IO^^Mq. This fraction is about 30% for filaments in the same mass range, and is 
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Fig. 3. — The upper left, upper right and bottom left are the comoving number densities 
(dashed lines) and mass fractions (sohd hues) for sheets filaments and halos respectively. 
The bottom right of the figure compares the comoving number densities of sheets (dotted 
lines), filaments (dashed lines) and halos (solid lines). The top sets of curves are for objects 
with masses of 10^*^ Mq, and subsequent curves show masses which are larger by 1 dex until 
the bottom sets of curves, which are for masses of IO^^Mq. 



only 13% for halos. At the present time, more than 99% of the cosmic mass is contained in 
sheets more massive than 10^°Mq, in contrast to viriahzed halos which, in this mass range, 
contain only about half of the cosmic mass. At 2; ~ 3, more than 10% of the cosmic mass 
was already assembled in sheets more massive than that of the Milky Way halo (~ IO-'^^Mq), 
while only 2% was assembled to halos in this mass range. Even at 2; ~ 10, more than 1% of 
the cosmic mass was already in sheets more massive than 10^*^ Mq, whereas the mass fraction 
in halos in the same mass range was ten times smaller. This suggests that the structure at 
high z is dominated by sheets. 
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At any given redshift, we may define a characteristic mass, M^{z), so that a{M^) is equal 
to the critical overdensity Ssc{z). In our ACDM model, the value of drops rapidly with 
redshift, from just over lO^^M© at z = to ~ IO^^-^Mq at z = 2. Objects at z whose 
mass satisfies a{M) = N6sc{z) = Na{M^,), where = 1,2, ■ ■ -, can be thought of as being 
increasingly unusual compared to objects of the characteristic mass M^,(z). The top two 
solid curves in the first three panels of Figure 3 show the abundance-redshift locii of sheets, 
filaments and halos whose masses satisfy a{M) = cr(M*) and cr(M) = 2a (M^). The panels 
show that the mass fraction in la sheets, 0.59, is larger than that in la filaments (0.32) 
and than that in la halos (0.17). A similar scaling applies to the more massive 2a objects, 
implying that, at any redshift, the typical structure is dominated by sheets. 

The bottom right panel in Figure 3 compares the number densities of sheets (dotted), fila- 
ments (dashed) and halos (solid) at different redshifts in the three other panels directly (i.e., 
these were the dashed curves in the other three panels): top-most curves show our model 
for masses of IO^^Mq, bottom-most for IO^^Mq, and the curves in between show results in 
which the mass changes by 1 dex. These curves show that, at masses larger than IO^^Mq, 
there are more sheets than filaments, and more filaments than halos, for almost all redshifts. 
However, at small masses and late times this trend begins to reverse: at z = there are 
fewer sheets than filaments with masses smaller than IO^^Mq, and there are more halos than 
sheets of mass IO^^Mq. Based on these mass functions, we can estimate the average mass 
of objects more massive than some given minimum mass. For instance, when the minimum 
mass is IO^^Mq then, at the present time, this average mass is 3.9 x IO^^Mq for sheets, 
2.1 X IO^^Mq for filaments, and 1.6 x lO^^M© for halos. The average mass of sheets is about 
2.5 times that of halos, again suggesting that large scale structure is dominated by sheets 
rather than by virialized halos. 



2.6. Conditional mass functions 

Halo abundances are expected to correlate with the overdensity of their surroundings — 
massive halos populate dense regions (Mo & White 1996). However, our model predicts 
that halo abundances will also correlate with the morphology of their surroundings. To 
illustrate, the solid curve in Figure 4 shows the mass fraction of IO^^Mq sheets that is in 
filaments of mass m at z = 0. The figure indicates that the most probable filament mass, 
7 X IO^^Mq, is a substantial fraction of that of its parent sheet. Although we do not show it, 
the analytic estimate of this quantity, given by inserting our barriers for filaments and sheets 
in equation (7) of Sheth & Tormen (2002), is in excellent agreement with this distribution. 
This solid curve should be compared with the dashed one, which shows the mass fraction 
in halos (rather than filaments) of mass m that are known to be in sheets of mass lO^^M© 
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Fig. 4. — Mass fraction of IO^^Mq sheets that is in filaments (sohd) and halos (dashed) of 
mass m all at z — 0. Dotted curve shows the mass fraction of IO^^Mq filaments at 2; = that 
is in halos. Dot-dashed curve shows the mass fraction contained in halos of mass m within 
an average volume of the universe of the same mass (IO^^Mq). The differences between the 
dotted and dashed curves indicate that, at fixed large-scale overdensity, the halo population 
is expected to be correlated with the morphology of the surrounding large-scale structure. 

at 2; = 0. (In this case, the analytic estimate is not as accurate.) Notice that the most 
probable halo mass is significantly smaller than that of the parent sheet, and the halo mass 
function in sheets is skewed significantly towards lower masses than that of filaments. As 
a final comparison, the dotted fine shows the mass fraction of 10^^ M© filaments (rather 
than sheets) that is in halos of mass m, and the dot-dashed curve shows the mass fraction 
contained in halos of mass m within an average volume of the universe of the same mass 
(10^'^ Mq). Here again, the mass function within sheets and filaments is biased compared to 
the average universe and the most probable halo mass is a substantial fraction of the mass 



12 



of the parent filament. Tfie significant differences between tfie dotted and dashed curves 
illustrate an important physical implication of our model: if the properties of a galaxy are 
correlated with the mass of its host halo, then the galaxy population in sheets of a given 
mass is expected to be different from that in filaments of the same mass. In other words, at 
fixed large-scale mass, galaxy properties are expected to be correlated with the morphology 
of the surrounding large scale structure. 
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Fig. 5. — Conditional distribution of sheets centred on a smaller scale filament of mass 
lO^^M© (solid). Dotted and dashed curves show the distribution of filaments and sheets 
which surround IO^^Mq halos. 

Figure 4 shows the distribution of masses of smaller scale objects when the large scale mass 
and morphology is fixed. But our model also allows one to estimate the mass and morphology 
of the environment which surrounds a given smaller scale structure. To illustrate. Figure 5 
shows conditional mass functions in which the small-scale object is fixed, and we study the 
distribution of the surrounding large scale mass and morphology. The solid curve shows 
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the distribution of sheets around IO^^Mq filaments, and the dotted and dashed curves show 
the distributions of filaments and sheets around lO^^M© halos. (Once again, the analytic 
estimates of these distributions are in reasonable but not perfect agreement.) Notice again 
that the mass of a sheet is well correlated with the mass of the filament it surroimds, as 
are the masses of filaments with those of the halos they surround. On the other hand, 
halo and sheet masses can be very different. These scalings indicate that our model can be 
used to estimate how the surrounding large scale structure correlates with different galaxy 
populations. 



Our model can be used to provide a simple estimate of how halos grow. To see why, consider 
a random walk with height equal to the rms fluctuation: a. The typical mass scales of sheets, 
filaments and halos are given by setting this value equal to the various barrier heights: a — 
(^eci,2,3(c, z)- To see what this implies, approximate the three moving barriers in equation (2) 
as 



Then the typical mass scale of filaments is (Jf ~ Sc{z). Similarly, a characteristic mass scale 
for sheets is CTg ~ Sc{z) — 0.5 CTg, implying cTg = (2/3) at, whereas the mass scale of halos is 

CTh = 2(Tf. 

Massive objects are associated with more extreme fluctuations. We can approximate by 
studying the case in which the walk height is twice the rms value (a 2a fluctuation). In this 
case, o-f = Sc{z)/2, = (2/5) 6c{z) = (4/5) cxf, and (Xh = (2/3) 6c{z) = (4/3) (Tf. Notice that, 
for massive objects, these factors are substantially closer to unity than they are for more 
typical objects. This suggests that the most massive halos are a substantial fraction of the 
mass of the filaments they populate, and these filaments are a substantial fraction of the 
sheets in which they are embedded. Thus, massive halos will appear to be accreting most of 
their mass from filaments. In contrast, a lower mass halo may be substantially less massive 
than the filament it inhabits. The growth of such a halo will appear to be dominated by 
mergers with the other small halos which populate the same filament (and possibly the other 
nearby filaments). In this respect, our model is qualitatively consistent with observations 
which suggest that massive clusters grow by accretion along filaments, whereas the growth 
of lower mass objects is less obviously anisotropic. 



2.7. Filamenteiry versus isotropic mass growth 




Sec2{cr, Z) 
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2.8. Identification of sheets, filaments and halos in the cosmic web 



Although our approach provides a framework for discussing the morphology of the cosmic 
web, we have, so far, not made any statements about precisely how the sheets, filaments and 
halos in our formalism are to be identified in cosmic density fields. Our collapse model is 
cahbrated so that the object which forms after collapse has been completed along all three 
axes, i.e., a halo, has the same density relative to the background as is expected in the 
spherical collapse model. Sheets in our model should be approximately this overdensity to 
the one-thirds power, whereas filaments arc this to the two thirds power. Although, formally, 
the density of a halo depends on the background cosmological model, it is common practice 
to identify halos at a given epoch as objects which are 200 times denser than the critical 
density at that epoch. Therefore, filaments and sheets should be approximately 36 and 6 
times denser than the critical density. This provides a simple rule of thumb for identifying 
sheets and filaments in the dark matter density field. Note that, in the context of our model, 
it may be more appropriate to think of massive halos as lying at the centers of filaments, 
rather than defining their endpoints. It will be interesting to see how these simple estimates 
compare with measurements in simulations. 



3. Discussion and Conclusion 

We extend the excursion set approach to quantify how the cosmic web is made up of sheets, 
filaments and halos. Our model assumes that objects form from a triaxial collapse; we define 
sheets as objects which have collapsed along only one axis, filaments as objects which have 
collapsed along two axes, and halos as objects where all three axes have collapsed. Therefore, 
our model requires specification of exactly how triaxial collapse occurs. Appendix A discusses 
our preferred collapse model, compares it with Zeldovich's approximation, and shows how the 
analytic arguments of White & Silk (1979) can be extended to provide an accurate analytic 
description of our ellipsoidal collapse model. 

The details of how a patch collapses depends on its initial overdensity 6, and on its initial 
shape parameters e and p. To embed this collapse model in the excursion set approach 
requires study of three-dimensional random walks crossing a barrier B{6,e,p). Here, we 
follow Sheth et al. (2001) and study the simpler problem in which e and p are replaced by 
representative values, and then study a one-dimensional boundary crossing problem. We 
emphasize that this is only an approximation, albeit a useful one. 

For any redshift z, insertion of the representative values of e and p in our collapse model 
provides estimates of the critical overdensities required for collapse along one, two and three 
axes. We find that these overdensities, Seci{z, a), Sec2iz, cr) and Secsiz, a) depend on both time 
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and mass (equation 2). Because of the dependence on cr, in the language of the excursion 
set approach, each of these critical overdensities is a 'moving barrier' (Figure 1). Insertion of 
each moving barrier into the excursion set approach provides estimates of the mass fraction 
in sheets, filaments and halos as a function of mass and time (Figure 2). 

In our model, halos of a given mass ruh populate filaments which are more massive, and the 
filaments themselves are surrounded by even more massive sheets. Hence, the characteristic 
masses of sheets are predicted to be substantially larger than of filaments or halos. Every 
halo at a given time was previously a filament of the same mass, and before that, a sheet. 
A halo of mass IO^^Mq today was a filament at redshift z ~ 0.45 and a sheet at 2; ~ 1 
(c.f. equation 2). Halo abundances are expected to correlate with the overdensity of their 
surroundings (massive halos populate dense regions). Our model predicts that, at fixed 
large-scale overdensity, halo abundances will also correlate with the morphology of their 
surroundings (Figures 4 and 5). Therefore, in models where the properties of a galaxy 
are correlated with the mass and formation history of its host halo, our model provides a 
framework for describing correlations between galaxy properties and the morphology of large 
scale structure. For instance, distributions like those in Figure 4 may be used to study if 
morphological structures in the galaxy distribution, like the SDSS Great Wall at 2; ~ 0.08, are 
unusual. And distributions like those in Figure 5 provide a framework for understanding if 
galaxies at high redshift form preferentially in sheets or filaments (Mo et al. 2005). Because 
our model exhibits various physically appealing features, we anticipate that it will provide a 
useful framework for quantifying the relation between galaxies, halos, and the cosmic web. 
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A. Appendix 

This Appendix describes the triaxial collapse model we use in the main text. Although 
the evolution must be solved numerically in general, we also discuss a reasonably accurate 
analytic approximation to the evolution, and show that it is considerably more accurate than 
the Zeldovich approximation. 



Let Ak denote the scale factors for the three principal axes of the ellipsoid. Then 



-AnGpAk 



(Bond & Myers 1996), where p oc is the mean density of the universe. 



(Al) 



p ^ ^ 

is the relative overdensity, and and Ag^tfe denote the interior and exterior tidal forces. 

In particular, 

dr 



2 

b'kit) ^bkit)--, where bkit) = Ai{t)A2{t)As{t) J 



(A3) 



This preprint was prepared with the A AS IM^jX macros v5.2. 
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and the linear approximation for the external tides is 



(A4) 



where = a(tj). Note that the A^s are the initial eigenvalues of the strain tensor; they are 
related to the initial density contrast 6i and the shear ellipticity e and prolateness p by: 

Ai = (5,/3)(l + 3e + p), A2 = ((5,/3)(l-2p), and A3 = (V3)(l - 3e (A5) 



The initial conditions are set by Zeldovich approximation: 

Akik) = ai[l - Xk{ti)], and Ak{ti) = H{ti)Ak{ti) - aiH{ti)\k{ti 



(A6) 



These expressions reduce to those of the spherical collapse model if the two tidal force terms 
are not included. Note that, as a result of the tidal forces, an initially spherical region will 
be distorted into a collapsing homogeneous ellipsoid. 

The tidal gravitational forces accentuate the asymmetry of the ellipsoid, so that the collapse 
of the three axes can happen at very different times. The shortest axis will collapse first, 
followed by the intermediate axis and then by the longest axis. Absent angular momentum, 
the collapse of any given axis continues to arbitrarily small sizes, a well-known feature of the 
spherical collapse model as well. 

Given initial values of 5j, ellipticity e, prolateness p, and epoch Oj, equation (Al) must be 
solved numerically for each axis A^. However, it is straightforward to extend the analytic 
approximation for ellipsoidal collapse provided by White & Silk (1979) so that it reduces 
self-consistently to the Zeldovich approximation at early times. In particular, we write 



Ak{t) 



a{t) 



Ak{ti) l-D{t)Xk 



a{t) 



Ah{ti) 



D{t)di ae{t) 
3 a{t) 



(A7) 



where Afi{ti) — 3/ ^^(ij)"^, D{t) is the linear theory growth mode, and ae{t) is the 
expansion factor of a universe with initial density contrast 5i = Afe(tj). Note that the first 
term in the expression above is the Zeldovich approximation to the evolution. Also note that, 
if all three axis are initially the same, then all the A^s are the same, and A = Si/3. In this case, 
the perturbation is a sphere, and the expression above reduces to A{t) —>■ A{ti) [ae{t)/ai], so 
the approximation (A7) is exact. 

Figure 6 compares the full numerical evolution of the three axes in this model (solid) with the 
Zeldovich approximation (dotted) and with equation (A7) (dashed), when the initial values 
are {S,e,p) = (2/400,0.2,0) in an Einstein de-Sitter universe. This shows clearly that, at 
late times, equation (A7) is significantly more accurate than the Zeldovich approximation. 



19 




50 100 150 200 250 300 

a/a. 



Fig. 6. — Comparison of full ellipsoidal collapse model (solid), Zeldovich approximation 
(dotted) and the approximation in equation (A7) for e = 0.2 and p = in an Einstein de- 
Sitter model. The three sets of curves show the comoving lengths of the shortest, intermediate 
and longest axes. These correspond to the formation of sheets, filaments and halos. 

The numerical solution allows each axis to shrink to arbitrarily small sizes. In studies of 
the spherical evolution model it is standard to assume that virial equilibrium leads to a 
non-negligible freeze-out radius. For an Einstein de-Sitter model, this freeze-out radius is 
Ak — 0.177a, so the final object is IStt^ 178 times denser than the background universe. 
Bond & Myers (1996) suggested that this same factor could be used in the ellipsoidal collapse 
model as well: hence, we stop the collapse of axis k by hand when A^/a — 0.17. 

Thus, in this model, initial values of and yield initial values of 5j, e and p. The collapse 
model then provides an estimate of the times at which each axis freezes-out. In particular, 
note that if e and p are given, then there is a unique value of 5i which will produce collapse 
of the kth axis at redshift z. The symbols in Figure 7 show the critical densities for collapse 
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Fig. 7. — Dependence of barrier height on mass for collapse along one, two and three axes 
(bottom to top). Symbols show the relations determined from the eUipsoidal collapse model, 
and curves show the analytic forms given in the main text (equation 2). 

at 2; = when p = and e = {a/6)/V5 along one (bottom), two (middle) and three (top) 
axes. In the random walk model discussed in the main text, these values of e and p, and this 
relation between the critical density required for collapse at z and the initial shape play a 
central role. The solid curves show the simple fits to these relations given in equation ( 2). 

Notice that the critical density for collapse of the second axis is very similar to that for a 
spherical model. The approximation (A7) provides an easy way to see why. First, note that 
to second order in 5j, 

M^~l_^_2(3e2+p2) 
3 

Second, recall that when p = then A2 = Si/3. Third, when p = and e = {a/6)/^/5 (our 
representative values) then the final term in the expression above is (2/15)aj^, so Afi{ti) ^ 
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A(ti) - Qi (2/15)al Therefore, for o", < 1, 



1 - D{t) A2 



a{t)Ah{ti) 
aiA2{ti) 



1 - 



3 a{t) 



(A9) 



in this approximation, the second axis evolves exactly as in a spherical model with initial 
over density S^. 
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